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We test a non equilibrium approach to study the behavior of a Bose-Fermi mixture of alkali atoms 
in the presence of a Feshbach resonance between bosons and fermions. To this end we derive the 
Hartree-Fock-Bogoliubov (HFB) equations of motion for, the interacting system. This approach 
has proven very successful in the study of resonant systems composed of Bose particles and Fermi 
particles. However, when applied to a Bose-Fermi mixture, the HFB theory fails to identify even 
the correct binding energy of molecules in the appropriate limit. Through a more rigorous analysis 
we are able to ascribe this difference to the peculiar role that bosonic depletion plays in the Bose- 
Fermi pair correlation, which is the mechanism through which molecules are formed. We therefore 
conclude that molecular formation in Bose-Fermi mixtures is driven by three point and higher order 
correlations in the gas, unlike any other resonant system studied in the context of ultra-cold atomic 
physics. 



I. INTRODUCTION 

Feshbach resonances have been recently discovered 
in ultracold mixtures of bosonic fermionic alkali atoms 
pi 0. Together with the achievement of degenerate 
states of such systems 0. [J. . this experimental feat has 
opened investigative opportunities for the study of new 
ultracold regimes. From the theoretical point of view, 
on the other hand, studies of Bose-Fermi mixtures to 
date have been mostly limited to non resonant physics, 
focusing mainly on mean field effects in trapped sys- 
tems Jail llEl IH IH 03, phases in optical lattices 
0> lloT ll6lTl7|. or equilibrium studies of homogeneous 
gases, focusing mainly on phonon induced sup erfluidity 
or beyond-mean-field effects [H El El El El 

This paper introduces a time dependent theory of the 
Bose-Fermi mixture that accounts for the resonant in- 
teraction. In systems where the resonant interaction is 
between two boso ns ]24l Eol E(j or between two fermions 
ElElElEllIlllJthe theory of "resonant superfluid- 
ity" has already been articulated. This theory is, so far, a 
big success. In the Bose case, it quantitatively describes 
the coherent conversion of bosonic atoms into bosonic 
molecules and back. Indeed, Ramsey interferometry on 
this system, coupled with this theoretical analysis, has 
produced the most accurate interaction potentials yet be- 
tween ultracold rubidium atoms |£j|. In the Fermi case, 
the theory has produced important qualitative insights 
into the crossover regime between weakly-interacting 
Cooper pairs on the one ha nd, and Bose-condensed 
molecules on the other [IlElElEIIllIlEllIil. 

It seems worthwhile, therefore, to adapt the same level 
of theory to the resonant Bose-Fermi mixture. In this 
paper we formulate the problem by writing down the rel- 
evant equations of motion at the level of Hartree-Fock- 
Bogoliubov (HFB) approximation. The equation of mo- 
tion are suitably number- and energy-conserving, as are 



their counterparts in boson of fermion systems. How- 
ever, in sharp contrast to these systems, the HFB theory 
applied to the Bose-Fermi resonance does not provide 
quantitatively reasonable results. Specifically we show, 
by direct numerical solution that the theory cannot re- 
produce the binding energy of a Bose-Fermi molecule, 
even in the limit of low density. 

The source of this difficulty lies in the approximate 
treatment of three-body correlations in the theory. The 
molecules, after all, are composed of two atoms, so 
the atom-atom-molecule correlation function is of cen- 
tral importance in determining properties of the result- 
ing molecules. In the HFB theory, this three-body cor- 
relation function is approximated in terms of two-body 
correlation functions, which is adequate for Bose-Bose 
and Fermi-Fermi resonances, but not for the Bose-Fermi 
mixture. Ultimately, the critical missing piece will turn 
out to involve the non condensed bosonic atoms. 

This paper is organized as follows: We begin our dis- 
cussion in section [H] by introducing the Hamiltonian of 
the system, and justifying such choice. We then proceed 
to outline the Bogoliubov-Born-Green-Kirkwood-Yvon 
(BBGKY) formalism used to derived the HFB equations 
of motion, and show the form they take in free space. In 
section ITTT1 we present our results, by first analyzing the 
equations by physical and analytical insight, and then 
presenting numerical results in support of our conclu- 
sions. 

Section ?? approaches the problem from an alterna- 
tive, perturbative point of view, relevant to low fermionic 
densities. From this analysis it is clear that molecular 
binding energies will not be recovered without adequately 
accounting for the non condensed bosons, thus pointing 
to their need for a higher-order theory. 
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II. THEORETICAL FORMALISM 



where 



A. The Hamiltonian 

We are interested primarily in the effects of resonant 
behavior on the otherwise reasonably understood prop- 
erties of the system. To this end we use a model which, 
in the last few years, has become one of the standards in 
the literature, and which was used to study the effect 
of resonant scattering in systems composed of bosons 
HHI3, and fermions Hi |H |H ^ IH H2 . Be- 
cause there is already a significant literature explaining 
the details involved in the choice of the appropriate model 
Hamiltonian, we only outline the extent of the approxi- 
mation involved in such a choice. 

An accurate approach to the problem would have to 
incorporate several scattering channels, since the reso- 
nance in question is a consequence of the intertwined 
behavior of the complex internal structures of the atoms. 
In a field theoretical sense, that would imply having to 
consider vector fields for the bosons and fermions with 
as many components as there are spin states involved 
in the interaction, and a non local interaction tensor of 
adequate size, to account for all coupling between such 
components. Fortunately, if we assume that the reso- 
nances in the system are sufficiently far from each other, 
such that it is possible to define a "background," or away 
from resonance behavior, we can focus on only one reso- 
nance at a time, which in turn makes it possible for an 
effective two channel model to describe the resonance. 
Furthermore, since the closed channel threshold is en- 
ergetically unaccessible at the temperatures of interest, 
we can "integrate out" the closed channel components of 
the fermion and boson field, in favor of a fermion fields 
which we identify as representative of the motion of one 
boson and one fermion, and which we dub the "molec- 
ular field." In the appropriate limit the molecular field 
identifies bound states between fermions and bosons. We 
emphasize that the molecular field is a theoretical arti- 
fice that alleviates the need to treat relative motion of 
two atoms on the natural scale of the interaction (tens 
of Bohr radii). However, this model is appropriate for 
the study of the systems at hand, typically composed of 
10 12 40 K atoms per cubic centimeter, whereby the char- 
acteristic length scales associated with the many body 
system are of the order of the inverse Fermi wavenumber, 
(thousands of Bohr radii), and the average interparticle 
spacing (tens of thousands Bohr radii), which is given by 
(ifp) 1 ^ 3 ' wnere P i s t ne atomic density. Lastly, since the 
coupling terms in the Hamiltonian represent an effective 
interaction, we can choose its functional form, and we do 
so by choosing to deal with contact interactions, which 
simplify the calculations immensely. 

The resulting Hamiltonian has the following form: 
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Here a p ,b p , are the annihilator operators for, respec- 
tively fermions and bosons, Cn is the annihilator operator 
for the molecular field [23, |3l], [33 ; 7 = iirab / mb is the 
interaction term for bosons, where a& is the boson-boson 
scattering length; and Vb g ,v, and g are parameters re- 
lated to the Bose Fermi interaction, yet to be determined. 
Also we define single particle energies e a — p 2 /2m a , 
where m a indicates the mass of bosons, fermions, or 
pairs, and V as the volume of a quantization box with 
periodic boundary conditions. 



B. Two Body Scattering Parameters 

The first step is to find the values for Vb g , u, <?, in terms 
of measurable parameters. We will, for this purpose, cal- 
culate the 2-body T-matrix resulting from the Hamilto- 
nian in eq. Integrating the molecular field out of the 
real time path integral, |4lJ leads to the following Bose- 
Fcrmi interaction Hamiltonian 



H 



2body 
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(3) 



This expression is represented in center of mass coordi- 
nates, and E is the collision energy of the system. From 
the above equation we read trivially the zero energy T- 
matrix in the saddle point approximation 



T = (V bg ~ 9 ~), 



(4) 



which corresponds to the Born approximation, and we 
proceed to match it to the conventional parameterization 
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H — H + Hi, 



(1) 



where ab g is the value of the scattering length far from 
resonance, As is the width, in magnetic field, of the res- 
onance, to;,/ is the reduced mass, and Bq is the field at 
which the resonance is centered. 

The identification of parameters between eqns, and 
|J5J proceeds as follows: far from resonance, \B — Bq\ >> 
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As, the interaction is defined by a background scatter- 
ing length, via Vb g — To relate magnetic field 
dependent quantity B — Bq to its energy dependent ana- 
log v, requires defining a parameter 5b = dv/dB, which 
may be thought of as a kind of magnetic moment for the 
molecules. It is worth noting that v does not represent 
the position of the resonance nor the binding energy of 
the molecules, and that, in general 5b is a field-dependent 
quantity, since the molecular binding energy approaches 
threshold quadratically. For current purposes we iden- 
tify 5b by its behavior far from resonance, where it is 
approximately constant. Careful calculations of scatter- 
ing properties using the model in eq. @, however, leads 
to the correct Breit-Wigner behavior of the 2-body T- 
matrix, as we show in section llVl 

Finally we get the following identifications: 



g = ^/V bg 5 B ^B 
u = 5b(B-B ). 

(6) 

For our calculations we use the 51 1G resonance in the 
40 K- 87 Rb system, the parameters we use in the calcula- 
tions to follow are ab g = — 202ao, 5b = 5.1 x 10~ 5 K/G, 
and A B = 1G. 



C. The Formalism 

We now move on to the many body analysis, and de- 
rive the Heisenberg equations of motion for the many 
body system. The way this is done, is to find equation of 
motion for correlation functions, f s (xi, ■■..x s ), which rep- 
resent the probability of finding s particles at positions 
X\,...,x s . As it turns out, the equation of motion for 
the correlation function fx will depend on the function 
/2, which in turn will depend on /s, and so on all the 
way to /at, where N is the total number of particles in 
the system. This is known as a Bogoliubov-Born-Green- 
Kirkwood-Yvon (BBGKY) hierarchy . In practice 
we will be concerned with momentum space correlation 
functions, but the idea is the same. 

Given the large number of particles in the system, it 
is impossible to calculate equation of motions for all cor- 
relation functions, and we need to invoke an approxima- 
tion. In practice, correlation functions are often calcu- 
lated only up to two-body correlations, s = 2. This is 
justified under the assumption that interactions are suit- 
ably "weak." Higher-order correlations are included in an 
approximate way by considering, not the actual atomic 
constituents, but rather combinations called quasiparti- 
cles. The quasiparticles are defined to be noninteract- 
ing, so that their higher-order correlation functions can 
be written in terms of second order correlation functions 



Using this qualitative idea we proceed to develop a 
more formal understanding. In statistical field theory, 
given an operator O, and Hamiltonian H, we define the 
thermal average of O with respect to H as < O >h= 
l/ZTr{Oe- 0H ), where /3 = (fesT)" 1 is the inverse tem- 
perature, and Z = Tr {e~@ H } is the partition function. 
In this framework, the 1-particle correlation function is 
defined as the thermal average of the number operator, 
with respect to the Hamiltonian of the system. 

In the quasi-particle representation, we define the anni- 
hilation operator for quasi-particles as a, reminding our- 
selves that it is a complicated function of a, a' , 6, w, c, c'. 
In momentum space, the 1-particle correlation function 
in this representation, will then be < a p \a P 2 >h , where 
H qp is the (noninteracting) quasi-particle Hamiltonian. 

Now we introduce the real approximation, namely that 
the quasi-particles can be written as linear combinations 
of all possible products of two operators (except for av- 
erages involving one fermionic and one bosonic operator, 
which are easily shown to vanish) . The procedure is then 
to find the Heisenberg equations of motion for these pairs 
of operators, and then averaging over the quasi-particle 
Hamiltonian 

i h^- t <0> Hqp =<[0,H}> Hqp , (7) 

which being Gaussian allows us to invoke Wick's theorem 
to decompose all higher order correlations in 1-particle 
correlations, thus truncating the BBGKY hierarchy. 



D. The Equations of Motion 

Before generating Heisenberg equations, we need to 
take a little care in the treatment of the Bose field, 
to properly treat the condensed part. To this end we 
perform the usual separation of mean-field and fluc- 
tuations of the Bose field, substituting b ( the zero- 
momentum component of the Bose gas) with a c-number 
<p =< bo >H qp , and identifying it with the condensate 
amplitude, while < b p ^Q >h = are the fluctuations. 
We insert these definitions in the Hamiltonian in eq. J5J), 
then proceed to calculate commutators. 

Since we wish to limit our analysis to a homogeneous 
gas, we note that the correlation functions fi(x,x') can 
be written in terms of a relative coordinate y = x — x' . 
Thus in momentum space f\ (p) is the probability to find 
a particle with momentum p in the gas, or in other words 
it is the momentum distribution of the system. 

Having taken all appropriate commutators, and ap- 
plied Wick's theorem, (for more details on the procedure 
see |36j , or Appendix ^ for the derivation of a sample 
equation.), we obtain the following self consistent set of 
equations of motion for the system: 
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(8.a) 
(8.b) 
(8.c) 
(8.d) 
(8.e) 
(8.f) 



(8.h) 



where t]~b{p) =< fcj^Q&p^o >H qp is the momentum distri- 
bution of non-condensed bosons, and ps — J ^jP 2 Vb(p) 
is the density of non-condensed bosons; k-b{p) =< 
bp^obp^o >H qp is the anomalous distribution of bosonic 
fluctuations, and As — J -§j?iP 2 kb{p) the anomalous 
density. Similarly r) Ft M(p) are the fermionic and molec- 
ular distributions, Pm, F the densities, and k_f,m(p), and 
Aj? m the anomalous molecular and fermionic distribu- 
tions and densities. Finally t]mf(p) =< c p a p >H qp and 
k>mf(p) —< Cpdp >h are the normal and anomalous 
distribution for molecule-fermion correlation, with the 
associated densities pm F and Amf- 



III. ANALYSIS AND RESULTS 



Equations I|8.al8.i|l describe the complete self- 
consistent set of HFB equations for the resonant BF mix- 
ture. Inspection of these equations, however, allows us 
to simplify the set quite dramatically, without sacrific- 



ing almost any of the physics thereby contained. First, 
we notice that the evolution of the anomalous fermionic 
densities km F {p), kf{p), and km(p) is entirely decoupled 
from the evolution of all other quantities, and can there- 
fore be considered separately. This implies that, since 
we are mainly interested in the evolution of the normal 
densities, we can eliminate without approximation all the 
anomalous ones. 

The next thing we notice is that the evolution of the 
normal and anomalous bosonic averages is completely in- 
dependent of the resonant interaction, and is controlled 
only by the background interactions between bosons 
and with fermions. For typical background interaction 
strengths, and cold enough temperatures, it is well estab- 
lished that quantum depletion is minor, and the system 
is well described at the Gross-Pitaevskii level of approx- 
imation. 

We can therefore write the following reduced set of 
equations: 





ih—<t> = {V bgPF + j\(t)\ 2 )<j) + gp* MF (9.a) 

d 

n dt VF ^ = ^sm((f)7j MF (p)) (9.b) 
d 

KglVMip) = 2g ^m{(f>-q MF {p)) (9.c) 
d 

ifr-gj_VMF(p) = [4 - 4' ~ u + H S M 2 ] vmf(p) - g4>* (vf (p) - mi?)) • ( 9 - d ) 



r 



Together with the prospect of simulating time dependent experiments, such a set of equations allow us to calculate 
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many characteristics of the system, which we could use 
to understand further physics or, more importantly at 
this stage, to test the theory against our knowledge of 
the system in various limits. 

A relevant quantity we can calculate to this end is the 
binding energy of the molecules. This can be done by an 
instantaneous jump of the detuning from large and posi- 
tive values, where we know the equilibrium distributions 
very well, to some other arbitrary value. The system thus 
perturbed oscillates at a specific characteristic frequency, 
which identifies as the (unique) pole of the HFB many 
body T-matrix of the system. For negative detunings, 
as shown below, this pole corresponds to the binding en- 
ergy of the molecules, dressed by the interactions in the 
system. 

Figure^shows a representative example of time evolu- 
tion of the condensate population (number conservation 
guarantees that all three populations oscillate with the 
same frequency) under the conditions described above. 
In this particular example, at time t=0 the detuning 
is suddenly shifted to — 5.110~ 6 K, corresponding to a 
magnetic field detuning of approximatively .1G. The re- 
sponse of the population shows an envelope function, in- 
dicated by the gray shaped area, that arises from nonlin- 
earities in the equations of motion. The inset shows that 
under this envelope is a well defined sinusoidal oscillation. 

The nearly monochromatic character of the response is 
made clear by Fourier transforming the time dependent 
population. The Fourier Transform shown in the second 
panel of Figgis strongly peaked at 5.4 x 10 _6 iG Simi- 
larly, the position of the peak in the frequency spectrum, 
for different final detunings, should map the molecular 
binding energy as a function of magnetic field. 

Figure |2 shows the results obtained by this method. 
This plot represents the binding energy of the molecules, 
dressed by the interactions in the system. This dressing is 
expected to be weaker for smaller densities of atoms and 
molecules. In this limit, we should thus recover the two 
body molecular binding energy, which can be calculated 
quite accurately from two body close coupling calcula- 
tions (solid line in fig |2J) ■ Instead we see that the pole 
behavior approaches the bare detuning (dashed line in 
fig 121, indicating that the renormalization of the binding 
energy obtained at the presented level of approximation 
is inadequate to correctly include the two-body physics. 
This behavior is in sharp contrast to the Bose-Bose res- 
onant interaction, where the correct binding energy is 
preserved at the HFB level [24| . This is also true for the 
Fermi-Fermi case |3^ . 

This discrepancy is due to the fact that the creation of 
molecules requires the formation of correlations between 
bosons and fermions, which, as shown in the following, 
cannot exist if the density matrix is assumed to be Gaus- 
sian. Specifically what is required is a more careful con- 
sideration of the noncondensed bosons 
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FIG. 1: The top panel represents the time evolution of the 
population of condensed atoms after detuning is suddenly 
shifted from infinitely positive to -5.1 K ( — .1G magnetic 
detuning). The bottom panel shows the absolute value of 
Fourier transform of said time evolution, the main peak in 
this graph represents the computed value of the binding en- 
ergy, which we see is about 5.410 -6 K. The system under 
consideration is composed of fermionic densities of 10 12 cm -3 , 
for a constant density ratio of five bosons per fermion. 



IV. THE IMPORTANCE OF BOSONIC 
DEPLETION 

The reason for the failure of the HFB theory is not 
immediately clear from the theory itself. To bring out 
the inadequacy of this theory in the dilute limit, we now 
recast the problem in an alternative perturbative form 
that can reproduce the correct behavior in the two body 
limit. This path integral approach will also lay bare the 
role of noncondensed bosons. 

What we will see in the upcoming analysis may be 
qualitatively understood in the following simple terms. 
A molecule in the gas can decay into a pair of "virtual" 
(i.e. non energy conserving) atoms, which can then meet 
again and reform the molecule. The incidence of these 
events modifies the behavior of the molecule, and an ap- 
propriate treatment of these virtual excitations is there- 
fore necessary to correctly include the two body prop- 
erties of the molecules in the many body theory. In 
particular, the molecules can decay forming a virtual 
non condensed boson, and the contribution of this set 
of events to the physics of the molecules turns out to be 
very important. An appropriate theory would therefore 
consider the coupling of the molecules to non condensed 
bosons explicitly, which implies that one has to include 
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FIG. 2: Plot representing the poles of the scattering t-matrix 
for the 511G 87 Rb - 40 K Feshbach resonance. The dotted 
line represents the "bare" molecular detuning as a function of 
field, as defined in the text. The solid black line is the correct 
binding energy of the molecular state, obtained by means of 
full close coupling calculations, while the gray solid lines are 
the eigenenergies obtained from equations I9.al9.dl for differ- 
ent atomic densities. From top to bottom on the right the 
grey lines refer to fermionic densities of 10 10 cm _3 ,10 11 cm~ 3 , 
and 10 12 cm -3 , for a constant density ratio of five bosons per 
fermion. We note that for lower and lower densities the calcu- 
lated binding energy incorrectly approaches the bare detuning 
instead of the correct two-body binding energy. 



3-body correlation functions differently depending on the 
quantum statistics of the constituents. For a Bose- 
Bose mixture, the correlation function is approximated 
(schematically) by 

(bbm) S3 (m) (b(-q)b(q)) + 2(b) (bm) . (10) 

The first term of the right of this expression allows ex- 
plicitly for virtual bosonic pairs of arbitrary momentum, 
provided that the molecular field (m) accounts for most 
of the molecules, which is assumed to be the case. Sim- 
ilarly, in a mixture of distinct fermions, the correlation 
function reads 

(/i/ 2 m) « (mXAGO/aC-?)), (11) 

and the same argument applies, since the molecules are 
bosons. 

For the Bose-Fermi mixture, on the other hand, the 
correlation function would be approximated by 

(bfm) « (b) (fm) + (/> (bm) + (m) (bf) . (12) 
The required virtual atom-atom pairs would arise from 
the third term on the right-hand-side of this expression. 
However, these molecules are fermions, which have no 
mean field, (m) = 0. The only surviving term is then the 
first one, which accounts only for condensed bosons, and 
somehow correlates the fermionic atoms to the fermionic 
molecules. This is only an indirect way to get the bosons 
and fermions correlated. 



in the equations of motion three point averages, such 
as < cta_ p+g / 2 fr(p+(}/2) >■ Since the HFB theory dis- 
regards three point averages, it only contains molecule- 
atom-atom couplings of the form < c> q a q > <fio, where 
molecules can only decay forming a condensed boson. 
It is straightforward to see that the HFB theory treats 



1. Two- Body Scattering 

The perturbative analysis begins by recasting the 
Hamiltonian in eq. (0 in terms of a 2-body action, in 
center of mass coordinates : 



S[il>, V f ,<f>, f , f , £ f ] = S B [0, 0t] + Sf ty, tf] + S M [£, C f ] + S c bP, ^ , 0, f , e, £ f ] , (13) 
where the field <p represents the bosons, tp the fermions, and £ the fermionic molecules, and where 



^-qVu ,p' TO ,p 



p p.p'q 
p 

Suit, ^ = J ^ + ef + „) C, P U P 

p 

Scty^A^,^] = J^Y. < P -q4, P ' +q VVp'^P +9 J t^E(4p^p^p + c - c )' ( 14 ) 

p>p'q pi 
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where Ticu is the frequency associated with the motion of As before we will then proceed to integrate out the 
the various fields. molecular degree of freedom |4l| to get : 
I 



pp' 



(15) 



where E is the collision energy between the fermions and the bosons. We then undergo the inverse transformation to 
obtain: 



du ^ ( g 2 



s'^e\ =-/|££K+- — ? — :) £U 



cj,pSuj,p 



P \ 2(m b +m t ) 



(16) 



Here £' represents the new effective (i.e. primed) 
molecules. The first line of fig. El shows the diagrams 
describing the resonant collisions between bosons and 
fermions. Here the continuous lines refer to fermions, 
the squiggles to bosons, and the broken lines to effec- 
tive molecules. Since we are looking for poles of the S- 
matrix, we can disregard the trivial fermion and boson 
propagators, and proceed, as outlined in fig. [3] to calcu- 
late the renormalized propagator for denoted as M, 
represented there as a heavy broken line. This object 
coincides with the T-matrix of the system, and shares 
its poles. Using the definition of the retarded molecular 
self energy S M given in fig. [3J and calling the molecular 
propagator M (again for £'), we get the following Dyson 
series: 



T = M = M - M T, M M - Mq E Mq E Mo 

Mn - MnE M M, 



(17) 



where T is the T-matrix for the collision, and which has 
formal solution 
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FIG. 3: Feynman diagrams representing the resonant collision 
of a fermion and a boson. Straight lines represent fermions, 
squiggles bosons, and dashed lines represent the effective com- 
posite fermions 



T = M 



Mq 1 + E M ' 



(18) 



These quantities take the explicit form 
M {E) = (v bg 

T, M (E) =i 



9 



E-v 
dto dp 



3/2 

V2n V n 2 ' 



(19) 



where m^f is the boson- fermion reduced mass, and A is physical nature of the contact interactions; we will dwell 
an ultraviolet momentum cutoff needed to hide the un- 
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more on that shortly. Finally inserting eq. 1191 into eq. 
1181 we obtain the following expression for the T-matrix: 



T(E) 



3/2 



+ t- 



E 



A 



(20) 



To show that this expression correctly represents the 
two body T-matrix for two channel resonant scattering, 
we will compare it to the results we know from standard 
theoretical treatments, 38], which teach us that: 



T(E) 



(S(E)-l). 



m bf ^2m b fE 
where S(E) is the S matrix given by 



(21) 



S(E) = e" 2l \^ 



mhfE ah 



E-e Q + iTy/E , 



(22) 



Here r v E is the width of the resonance, eo is a shift asso- 
ciated with the detuning with respect to threshold of the 
resonance, and a b f is the s-wave scattering length for the 
boson-fermion collision; all of these quantities can be ex- 
tracted from experimental obscrvables, through accurate 
two-body scattering calculations. 

From the parametrization of the zero energy T-matrix 
in eq.|0, and the E — > limit of 121|) . we easily derive 
r = y / 2mBFa-BF^ti^B. With these definitions we can 
relate equations H21J) and J2DJ, to find a regularization 
scheme for the theory, by substituting the non observable 
parameters g, v, and V bg by the A dependent (renormal- 
ized) quantities g, D, and V bg , such that the observable 
T-matrix will not be itself A dependent 

Following j23| we compare equations (|2Tj|) and l(2*T|) . in 
the limit E —> 0, where we have (once we include the 
definitions of the bare quantities) 



Vbg 



1 



m b fA 



V bg £ 



(23) 



Since we have one equation and three unknowns, we will 
have to insert some physics in the system, analyzing it 
one limit at a time. The first limit is far from resonance, 
where v — > 00 



V,. 



1 



mbfAVbg 
=3 



(24) 



We are now left with the task of defining the resonant 
quantities, and we have no more leeway to make physi- 
cally motivated simplifications. The equations which re- 
main are ambiguous, which leaves us with a set of possi- 
bilities for the choice of g and v. One way is to proceed 
as follows: insert eq. I|24|) into l|23|) . and solve for D, to 
get 



3 2 (1 



rubfAVbg. /mbfA 



+ -)• 



(25) 



7T t: g 

From inspecting the above equation we can choose a 
definition of g, which will also imply one for v, and we 
get (reporting also ea. l|24|l for completeness) 



Vbg = V bg 



1 



1 - 



m bf AV bg 



9 = 9 



-. m b fAVb g 

r — 5 — 

7T 



v = v + gg 



m b fkV bg 



(26) 



Using these definitions of V bg , g and D, together with the 
policy of imposing A as the upper limit of momentum 
integrals, will guarantee that obscrvables will not depend 
on the choice of A, as long as it is chosen to be bigger 
than momentum scales relevant to experiment. 



2. Many-Body Generalization 

Generalizing the above treatment from two to many 
particles, we must now account for the fact that, in a 
many-body system the molecular self energy is modified 
by the environment. Unlike in the scattering problem, 
the procedure outlined in the previous section is only 
an approximation to the full many body problem, but 
as some evidence seems to suggest, a pretty good one 

MM. 



To perform this generalization, one needs to calculate 
the many body self energy using the many body free 
green functions D° MB (uj,p) and Gu B , respectively for 
bosons and fermions, defined as [4l| 



D° MB (oj,p) = \cb\ 2 + 



Tiuj — (Jp + + i0 + sign(e£ - /is) 



(27) 
(28) 
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where hb(f) is the bosonic (fermionic) chemical potential. The molecular self energy becomes 

X M B(n,p) = -ig 2 f ^^£>o(q,a,)Go(p-q,Q-u,) = 
J Zir (Zirr 



(2tt) 2 

+9 ip-kf ( 27r ) 3 ^ - £ q - £ p-q + + MB' ( ' 29 ' > 



the two terms in this expression represent contribu- 
tions from condensed (£g /B (£!,p)) and noncondensed 
(S^ /B (17,p))bosons, respectively. 

The (approximate) many body self energy in eq. 1291 
can be easily shown to reduce to its correct two body 
counterpart defined in eq. I|19|) . when the densities and 
chemical potentials are set to zero. However, if the contri- 
bution due to the non-condensed bosons is omitted, 
then S MB clearly vanishes in the two-body limit, <fr — > 0. 
There would then be no renormalization of the molecular 
propagator, and the pole of the T-matrix would coincide 
with the bare detuning, as shown in fig. [2J 

We remind the reader that neglecting the non- 
condensed component of the bosonic field was a perfectly 
well justified approximation of eas. l|8.a| l- H8.i|l . which im- 
plies that those equations are already inadequate to re- 
produce the two body binding energies in the low density 
limit. Indeed the part of resonant term in the Hamilto- 
nian containing the bosonic fluctuations vanishes accord- 
ing to Wicks theorem, since it is an average of a three op- 
erator correlation with respect to a density matrix which 
is Gaussian, in the HFB approximation. To correct this 
problem one should extend the HFB approximation and 
explicitly include three, and possibly higher, particle cu- 
mulants, finding some other way to truncate the BBGKY 
hierarchy. The subtleties involved in such a calculation, 
however, are many, and non trivial, and will be the sub- 
ject of further work. 



V. CONCLUSIONS 

We have performed a study the non-equilibrium be- 
havior in Bose-Fermi mixtures subject to an interspecies 
Feshbach resonance, Using the HFB approximation. We 
have found that this approximation is not adequate to 
describe the system, which is quite remarkable since it 
has become one of the standard approaches to resonant 
cold atom physics due to its successes in Bose gases and 
two component Fermi gases. 

The reason of this failure is found in the way in which 
the theory treats non-condensed bosons. This problem 
could be corrected by the explicit inclusion of three (and 
possibly higher) point cumulants, which will allow for 
a mechanism through which bosons and fermions could 
correlate to form molecules. This task, however is beyond 
the scope of the current investigation. 
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APPENDIX A 

In this appendix we will present a sample derivation of 
one of the equations of motion, namely that for < 5^6 >. 
Starting with the Hamiltonian in coordinate space 

H = J dx ^(x)T F (x)i(>(x) + 
/^.(,)T«(,W,) + 
J dx ^(x)T M (x)^(x) + 
dx\<f>(x)\* + 

V bg jdX \<f>( X )\ 2 Mx)\ 2 + 

g J dx (x)<p(x)<p(x) +c.c) 

(1) 

, where T a (x) is the kinetic energy of molecules, bosons 
or fermions. 

We then write the bosonic field in terms of its average 
and fluctuations around it 4>{x) — <fio(x) + S(x), where 
<fro is a complex number. Inserting this expression in the 
Hamiltonian, we get the following 

H =E Q + J dx^{x){T F {x) + V bg \^{x)\ 2 )^{x) + 
J dx 8\x)T B \x)5{x) + J dx C f (x)T M (x)£(x) + 
7 / dx (4\M^)\ 2 \S(x)\ 2 + <P*o(x) 2 6(x)6(x) + 

M^) 2 ^(x)s^x)) + 

J dx (cj)* (x)5(x) + (t> Q {x)S f (x)) x 

(i\ M x)\ 2 + v bg mx)\ 2 ) + 
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7 / dx (4>l(x)5^(x)5(x)5(x) +c.c) 



(2) 



dx S\x)5' i (x)6(x)S(x) 



V bg / dx \6(x)\ 2 \i>(x)\ 2 + 



g J dx [f(a:)t (0o(x) + 5(a:)) ^(s) + c. 



where i?o is a constant which depends on </>q, and its 
relevant for its motion, but does not contribute to that 
of flttf. 

The next step is to calculate the commutator 
[tft(z)(5(z'),-H], and to take its average, thereby obtaining 



<[8\z)8(z'),H] >= 



(T B {z')-T B (z)) <5\z)5{z') > + 

7 [2\Mz')\ 2 < SHz)S(z') > +4>l{z) < 5\z)5\z r ) > - 

2|0 o (z)| 2 < SHz)S(z') > ~^ 2 {z') < S(z)6(z>) >] + 

M*') Or > \Mz')\ 2 + v bg Jdx <5\z')^{x)i>{x) - 

^o(z) (7 < S(z) > \M?)\ 2 + V bg J dx < 5(z)^(z)i>(z) >j + 

7 J dx [^S(x)(< ^ (z)^')^ (a:)<y(a;)*(x) > - < 6\x)5(x)S(x)S^ (z)6(z') >) 

<M^)(< (5 t (z)5(z')(5 t (a;)(5 t (a:) ( 5(2;) > - < tf(x)rf [x)6{x)tf {z)S(z') >)] + 
7 (< <5 t (z) ( 5 t (z')<5(z')<5(z') > - < 8\z)5\z)8{z)5{z') >) + 
H fl (< S^z)S(z')^(z')^z') >-< 5\z)5{z')^{z)^z) >) + 

g (< COO^V)^*) > - < ^M^CO >) • 

I 



(3) 



The next step is to apply Wick's theorem to correlation all correlation functions of odd order will vanish. We then 
functions of three or more operators. This implies that get 



< [6\z)6(z'),H] >-- 



(T B {z')~T B (z)) <5\z)S{z')> + 
1 [2\U^')\ 2 <^{z)5{z')> + 

<f> 2 (z) < 5\z)5\z') > -2|0 o (z)| 2 < 5\z)8{z') > -0SV) < S(z)S(z') > 
7 (< 5\z)5\z') >< S(z)S(z') > +2 < 5\z')8{z') >< S^z)S(z') > - 
2 < 5(z)S^(z) >< S\z)5(z') >-< 5 t (z)«5 t (z) >< 5(z)6(z') >) 
Vbg (< tf{z)5{z') >< ^(z')^z') >-< SHz)S{z') >< cbHz)^z) >) + 



(4) 



In free space, <po becomes a constant, and all two point of z-z', so that in momentum space they become functions 
correlations, which are functions of z,z', become functions of a single momentum. We thus obtain eq. I|8.b|) . 
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